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RESEARCH MEMORANDUM 


LOCAL ISOTROPY IN TURBULENT SHEAR FLOW 
By Stanley Corrsin 

SUMMARY 


The mean strain rate in turbulent shear flow must tend to make the 
structure anisotropic in all parts of the spectrum. It Is argued here, 
however, that, if the spectral energy transfer process destroys orienta- 
tion, the Kdlmogoroff notion of local Isotropy can still be justified 
In spectral regions -where the local transfer time is shorter than the 
characteristic time of the gross shear strain. 


INTRODUCTION 


Recent measurements on the anisotropy produced by homogeneous strain 
of the (approximately) isotropic turbulence downstream of a grid (refs. 1 
and 2) suggest reexamination of Kolmogorof f ' s postulate of local isotropy 
In shear flow (ref. 3)- This postulate has been confirmed experimentally 
in the sense that isotropic-type behavior has been measured for the small 
"eddies" (refs. 4, 5* 6, and 7)- Up to the present time there seems to 
be no contradictory evidence. 

Nevertheless, as pointed out by Uberoi In reference 2, the ubiquitous 
strain rate of turbulent shear flow must act upon eddies of all sizes, 
tending to make the turbulence anisotropic at all wave numbers . The 
largest eddies must be nonisotropic in any case since their dimensions 
are comparable with the width of the mean shear zone. 

Kolmogoroff ' s idea apparently was that, in the nonlinear transfer 
of energy from small to large wave numbers, the necessary direction 
preference of the large structure gets lost. This seems quite plausible, 
but at least three other effects are involved: (l) The strain- induced 

tendency to anisotropy at all wave numbers, (2) the general tendency 
toward isotropy (equipartition?) evident In the absence of strain rate 
(see refs. 1 and 2) (this must also occur at all wave numbers), and 
(3) the viscous dissipation to internal energy, especially at very large 
wave numbers . 
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It is questioned whether the characteristic t ime s for (a) inertial 
transfer to’ higher wave number t , (b) component transfer at the same 

wave number t^, and (c) viscous dissipation at each step of the spectral 
energy cascade t c - are sufficiently short to forestall the mean strain- 
induced anisotropy. In inertial spectral domains if either r & or 

is short enough, local isotropy remains a plausible notionj in the domi- 
nantly viscous domain, t c is the critical time. The appropriate charac- 
teristic time for comparison is, of course, the inverse of the mean strain 


rate , or 



n\-l 


in boundary-layer-type flows with principal mean 


velocity field f}(x,y) . 


This work was carried out at the Johns Hopkins University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. The author would like to thank Mr. H. M. 
Fitzpatrick and Dr. C. C.'Lin for their helpful comments. 


SYMBOLS 


D 

E(k) 

k 



q. 


any characteristic length 
"three-dimensional" energy spectrum 
magnitude of wave-number vector 

Koltnogoroff wave number (inverse of Kolmogoroff microscale) 
characterizing viscous part of spectrum 

component wave number corresponding to Kolmogoroff wave-number 
magnitude 

magnitude of turbulent velocity, / u.2 + 


q' = ft? + v^ + ^ 

R Reynolds number 

■1 Q f 

R-x turbulence Reynolds number, -=== — 

A (S v 


Ui velocity vector 


U mean velocity in x-direction 

u,v,w turbulent velocity components in x-, y-, and z-directions 
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V any characteristic mean velocity component 

x,y,z Cartesian coordinate axes 

7 Heisenberg spectral constant 

A dissipation scale (Taylor microscale) 

v kinematic viscosity 

T a (k) spectral inertial transfer time 

•^(k) spectral time for approach to isotropy in absence of gross 

strain rate 

T c (k) viscous decay time 

$ rate of dissipation of turbulent kinetic energy per unit mass 

Subscripts : 

I,j,k vector directions 

CHARACTERISTIC SPECTRAL TIMES 


For a rough estimate of characteristic spectral times a discrete 
energy cascade in the manner of On sager (ref. 8 ) is visualized. The 
assumption of similarity in successive steps requires a geometric pro- 
gression} that is, at each step Ak k. Assume next that at each spec- 
tral jump In the cascade process all statistical orientation is lost} 
that is, the energy arrives at wave number k in isotropic condition, 
no matter how anisotropic it had become in the previous stage. Possibly 
by direct dimensional reasoning, Onsager defines the characteristic 
(inertial) transfer time per stage as 


T a (k) 



( 1 ) 


where k is wave-number magnitude and E(k) is a three-dimensional 
energy spectrum 


j(u 2 + v 2 + w 2 ) = | q 2 



E(k) dk 


( 2 ) 
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Here u, v, and w are orthogonal velocity fluctuation components. 
Physically this form can be deduced from 


T _ Kinetic energy per stage 
a Energy transfer rate 


(3) 


From equation (2) (and Ak » k) the numerator is obviously kE(k) . By 
analogy with the form of the product expressing the rate of transfer of 

Su t 

energy from mean shear flow to turbulence, = in the general 

Sx k 

turbulent energy equation (see ref. 9) , the denominator is expressed as 


(Spectral velocity)^ 
Spectral length 


( kE )^/ 2 

l/h 




Hence, following equation (l). 


T a( k ) 


kE 

l/k^E^ 


1 


Within each step of the cascade, the characteristic time t-^ for 

adjustment of an induced energy inequality among the three components 
is also an inertial phenomenon expressible in terms of pres sure -velocity 
correlations, hence for want of more detailed understanding it is 
as sinned that 


T b T a 


( 4 ) 


The viscous decay time . T c (k) of the energy at any stage is 
physically 


_ Kinetic energy per stage 


(5) 


G Viscous dissipation rate 
The denominator here, by analogy to the form of the dissipation rate in 


, aui/dii, 3uA , 

general v - — It— ^ - + — — *■ I, is 
ax j\ Sx j Sx i / 


(Spectral velocity) 2 
(Specral length) 


kE 


(i/k)‘ 


vk^E 


v 


v 
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where v is the kinematic viscosity. Hence 


T c (k) 


vk^E 


that is, 


T (k) = 1 (6) 

c vk^ 

Before returning to the question of local isotropy, it is instruc- 
tive to write out the two inequalities which identify the inertial and 
the viscous regions of the spectrum. The former region is temporally 
characterized by the fact that the energy in any stage jumps to the next 
before there has been time for appreciable viscous dissipation: 

r a (k) « T c (k) (T) 

The latter region consists of the cascade stages in which the energy is 
dissipated before there is time for inertial transfer: 


T a ( k ) >> T c( k ) 

Substituting equations (l) and (6) into these. 


Inertial region: 
Viscous region: 


if 

if- 


00 


» 1 


« 1 


This serves as a check on the r definitions since 
the spectral Reynolds number . 


if 


(Ml 


( 8 ) 


(9) 


is simply 


INERTIAL REGION 


With the assumed model a necessary condition for local isotropy In 
the purely Inertial spectral range Is that T a be much smaller than the 
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characteristic mean strain time, most simply Qj . In boundary- 

layer type shear flows (this Includes free flows in which Prandtl’s 

boundary-layer approximation is valid) ^ ^2 is a good approximation 

^ oy 

to the principal mean strain rate where t?(x,y) is the mean component 
of mean velocity and y is the Cartesian coordinate with Its highest 
gradient. Thus, inertial local isotropy can he expected at wave numbers 
for which 


1 « — — - ( 10 ) 

l/k?E(k) i |5 
^ oy 

In the inertial region a spectral behavior roughly like 

E(k) ~ (11) 


can be expected since this Is a formal result when local isotropy is 
presupposed (see refs. 3 and 8 ). Here 0 is the rate of dissipation 
of turbulent kinetic energy per unit mass . With this spectral form, 
r (k) decreases monotonlcally : 

cL 


T a (k) « 0" 1 / 5 k" 2 / 5 (12) 

If this expression were valid out to indefinitely large wave numbers, 
inertial local isotropy could he found far enough out the spectrum no 
matter how high the gross strain rate. There would always be a range 
of k for which 


r a (k) « 



( 13 ) 


Of course, equation (12) cannot apply" to indefinitely large wave 
numbers, so T a (k) will actually increase (eq. ( 28 ) or (30)). In any 
case, the inviscid form, equation (12), converts equation (l3) to 


t 2/3 » 1 0-1/ 3 g 

2 oy 


( 14 ) 
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In a real fluid, the inertial range can exist only at wave numbers 
much smaller than the inverse of the Kolmogoroff microscale, that is. 


k « 



(15) 


so that equations (ll) and ( 15 ) make a pair of necessary conditions for 
inertial local Isotropy. They can now be combined into a very crude 
Reynolds number criterion. 


Introduce the following symbol conventions: 
then a »> Pj (b) If a » 0 and p » 7 then 
suppose that a, » p implies a = Q(20P): then 



(a) If cc , 2 / 5 » P 2 / 5 
a, »» 7 . For example, 
e »» 6 implies 


With this representation. 


equation ( 1*0 can be written as 


k >» 



( 16 ) 


which can now be combined with equation ( 15 ) as 




(17) 


For a turbulent shear flow In which the total production rate for tur- 
bulent energy from mean flow energy Is of the same order as the total 
dissipation rate (they are exactly equal in pipe flow). 



But, empirically. 



according to measurements In a variety of shear flows, so 
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Furthermore, approximate 0 by its Isotropic form 


0 * 5v (21) 

A 2 

where A Is the Taylor microscale. Equations (20) and (21) convert 
equation (l7) to 



»»> 16 


( 22 ) 


In turbulent shear flows it seems reasonable to define the turbulence 
Reynolds number R^ as 


where 


x (5 v 



2 = q 2 = u 2 + v^ + w^ 
VISCOUS REGION 


(25) 


In the predominantly viscous region of the spectrum, 


k » 



iA 


(24) 


and the temporal inequality for possible isotropy is 


T c ( k ) 



1 

1 dU 

2 dy 


(25) 


or 




1_ dU 
2v 8y 


(26) 
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Since T c (k) decreases monotonically with k, it is clear that 
there will always be some wave number above which this necessary inequal- 
ity will be satisfied. 

It should be noted that in this spectral range the simple dimen- 
sional spectral transfer theory of Heisenberg in reference 10 gives 

E(k) « k“^ (27) 

a result which has had rough experimental confirmation. This leads to 
an increasing Inertial time 

T a (k) « k 2 (28) 

Although T a is of negligible dynamic significance in the viscous region, 
this suggests a detailed look at the mixed region, where 



Here both and t q are important and conceivably may violate equa- 

tion ( 13 ). 


MIXED REGION 


In the spectral region with both inertial transfer and viscous dis- 
sipation, -r (k) can be estimated from equation (l) by inserting for 

a 

E(k) the function obtained by solving Heisenberg's equation (see ref. 10) 


E(k) 




\9rj 

37 2 0 


-1-V3 


-5/3 


Here r c (k), as given by equation ( 6 ), is independent of E(k) . 


(29) 


The value of Heisenberg 1 s constant 7 in equation (29) has been 
variously estimated from experiment in the range 0.2 to 0 . 85 * with 0.4-5 
an acceptable compromise value. For simplicity, take 7 = 4/9. Assuming 
the "Kolmogoroff wave n umb er is written as 
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equation ( 29 ) turns equation (l) into 


/x\ - e-^V 17 ^ 1 / 2 


a Ik, 




i + 



(50) 


For values of k/k c sufficiently small that the second term in the 
brackets is negligible, this reduces to the inertial estimate, equa- 
tion (12) . 


The viscous decay time in terms of k/k c is 


= ^ 1/2 ( t )' 2 (31) 

Figure 1 is a dimensionless plot of T a and t . The values of 
r & and t c are necessarily of the same order in the vicinity of k » k c , 

since this is the region of equal inertial and viscous effects, that is, 
of Reynolds number about unity. 

By equating the derivative of equation (JO) to zero, the minimum 
value of r a is found to be 



T a (0.40) 



(52) 


On the whole, a necessary condition for local isotropy is that 


1.8 



« 


l 

1 £2 

2 3y 


(55) 


or simply 



(3fc) 


A less certain but more convenient criterion follows from the use 
of equations (20) and (21): 


R* » 15 


(55) 
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It Is worth, noting that for ky^, = 2, for example, the inertial 
o 

condition t_ « — — could he violated while the viscous condi - 

du/dy 

tion t_ « — - — is satisfied. Yet this is a spectral range in which 

inertial forces are not negligible, and equation (53) does not rule out 
the possibility of an anisotropic local spectral range for k> k„. Of 

course, for k » k c , t q is not pertinent and t c « (T a ) Tn ^ n , so equa- 
tion (55) implies isotropy here. 


WAVE NUMBER OF TURBULENT ENERGY PRODUCTION 


Following Prandtl's simple kinetic -theory type of model in turbulent 
shear flow, it is inferred that the production of turbulent energy from 
mean flow energy is primarily due to the lateral fluctuating motion of 
"fluid balls" in the presence of a velocity gradient. Therefore the mean 
production rate must depend at least upon the lateral root-mean- square 
fluctuation v* and the chief mean velocity gradient bu/by in a 
boundary-layer-type flow. 


Of course, the turbulent energy equation gives the rate per uni t mass 
as uiv Presumably this energy Is fed into the turbulence over a wide 

ay 

range of wave numbers, each receiving an allotment proportional to its 
spectral contribution to uv. 


It Is noted, however, that v* and bu/by are sufficient to repro- 
duce the dimensions of wave number k, so an estimate of the order of 
magnitude of the energy production wave number In, Is 



( 36 ) 


Clearly local Isotropy can be expected only at wave n umb ers much 
larger than kp : 


k» 4 


au 

by 


( 37 ) 


For turbulent shear flows whose gross production rate is of the same 
order as the dissipation rate, equations (20) and (21) can be Introduced 
so that 
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kX » 


30 

ST 


(38) 


For inertial local isotropy, l/k must also be much larger than 
the Kolmogoroff microscale; that is. 


k « 


iA 


which may be written together with equation (37) as 


or 


£ «« ft) 


iA 


(39) 


(40) 


(41) 


™ ^ 

Using equations (20) and (21) this gives a Reynolds number inequality 


3/2 


which is much like equation (22) . 


»» 25 


(42) 


COMPARISON WITH EXPERIMENT 


The three semiempirical Reynolds number criteria can be summarized 
as follows: 

(a) Lower bounds on R^ for the existence of an Isotropic inertial 
subrange are given by equation (22) based on Inertial transfer time and 
Kolmogoroff wave number and by equation (42) based on turbulent produc- 
tion wave number and Kolmogoroff wave number. 

(b) The lower bound on R^ for local Isotropy in. Inertial or mixed 

range is given by equation (35) based on minimum inertial transfer time 
assuming a Heisenberg spectrum. 

Next these inequalities are compared with the conditions in typical 
past experiments. In the round-jet experiments which include those of 

references 5 and 11 R^ «* 120, hence R-^/ 2 « 1,300. Since there is no 
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a priori basis for assigning a specific numerical factor to the strong 
inequality a » p, the experimental results are used. In these Jet meas- 
urements the shear correlation spectrum reached zero at a one-dimensional 
wave' number of about 45 per centimeter, whereas the Kolmogoroff wave num- 
ber k c (the inverse of the Kolmogoroff microscale) is 100 per cent ime ter. 

Since kc is a three-dimensional wave -number magnitude, comparison 
is properly made by defining a corresponding one-dimensional parameter k^. 

A plausible, though arbitrary, definition Is the first moment of the one- 
dimensional spectrum that corresponds to a Dirac-function three-dimensional 
spectrum, that Is, a spherical shell of radius k Q : 



Therefore local Isotropy is reached Just around the Kolmogoroff 
region where the spectral Reynolds number is of order unity. Hence 
equation (35) should be Just barely satisfied, and it is noted that 
this is so if a » p is interpreted to mean a 8p . 

In this flow there certainly can be no isotropic inertial subrange, 
and, in fact, neither equation (22) nor (42) Is satisfied. The largest 
value of R^ attained in a well-defined turbulent flow appears to be 

that in Laufer's pipe (see ref. 12), R^ «= 250. In this case, 

R-^/ 2 « 4,000, and he does find some evidence for an isotropic inertial 
subrange . 

In conclusion. It seems worthwhile to set up a turbulent flow with 
still higher values of R^ than those attained by Laufer, perhaps greater 

than -or equal to 500. Since R^ increases more slowly than R (where 
R = VD/v and V and D are characteristic gross velocity and width)’'' 
for fixed geometrical boundaries, values of 500 or higher will not be 
easy to attain with air In a small laboratory. 


The Johns Hopkins University, 

Baltimore, Md., August 1 6 , 1956 • 

■'"Neglecting the slow decrease in turbulence level which sometimes 

l/2 

accompanies increasing Reynolds numbers in shear flows, R^ cc R ‘ . 
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Figure 1.- Characteristic spectral times (assuming a Heisenberg 

spectrum) . 
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